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ABSTRACT 

A lifting surface theory has been developed 
for a helicopter rotor in forward flight for 
compressible and incompressible flow. The 
method utilizes the concept of the linearized 
acceleration potential and makes use of the 
vortex lattice procedure. Calculations demon- 
strating the application of the method are given 
in terms of the lift distribution on a single 
rotor, a two-bladed rotor, and a rotor with 
swept-forward and swept-back tips. In addition, 
the lift on a rotor which is vibrating in a 
pitching mode at 4/rev is given. Compressibi- 
lity effects and interference effects for a 
two-bladed rotor are discussed. 


INTRODUCTION 

Rotating lifting surfaces are an integral 
part of the propulsive unit of every aeron- 
autical and nautical vehicle, from the 
compressor and turbine blades of jet engines, 
the pumps for rocket engines, to propeller and 
helicopter rotors. The aerodynamics of these 
rotating elements has been under extensive study 
since the advent of the airplane and with a 
combination of experimental and analytical 
approaches, succcessful designs have been 
achieved. In many cases, two-dimensional theory 
has been used, usually modified by an assumed 
spanwise distribution, and inflow velocities. 
This paper presents a compressible, lifting 
surface method for a helicopter rotor in forward 
flight within the limits of linearized theory. 

The method is based on the concept of the 
acceleration potential, originally introduced by 
Kussner (1941). The method was first applied to 
an oscillating wing in uniform translatory 
motion including effects of compressible flow by 
Runyan and Woolston (1957). The acceleration 
potential approach has now become standard for 
the determination of the unsteady aerodynamic 
forces for flutter studies of lifting surfaces 
in rectilinear motion. 

The first use of the acceleration potential 
approach for a rotating system was made in a 
paper by Hanaoka (1962) for the loading on a 
marine propeller in incompressible flow. The 
acceleration potential has been used in the past 
in studying the propeller noise problem, but in 
all of these noise propagation cases the problem 
was specialized early in the analytical develop- 
ment to the so-called far-field case usually 


with a stationary observer, whereas the lifting 
surface theory is essentially concerned with the 
details of the near-field case for a co-moving 
observer as well as the satisfaction of certain 
of certain boundary conditions. Runyan (1973) 
utilized the acceleration potential approach to 
obtain a solution to the oscillating propeller 
in compressible flow. Dat, (1973), has derived 
a general expression for an acceleration doublet 
for any motion. Pierce and Vaidyanathan (1983) 
have treated the helicopter rotor in forward 
flight using the method of matched asymptotic 
expansion for the incompressible case. The 
procedure developed here involves the precise 
numerical integration over the surface of the 
rotor in a time frame. The method sets forth a 
formulation of a fundamental three dimensional, 
compressible, unsteady aerodynamic theory for 
propellers and helicopter rotors. 

The next section contains a brief 
derivation of the fundamental equations, includ- 
ing a discussion of some implications of the 
equations. The third section contains a 
description of the method of solution. Finally, 
the results of some calculations for the several 
specific examples are given. 
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SYMBOLS 

rotor blade area 
aerodynamic influence 
coefficients 
Fourier coefficients 
speed of sound 
chord of rotor 

thrust coefficient per blade 

(thrust/irpn 2 Rt 4 ) 

vector distance from doublet to 

downwash point 

absolute value of 5 

unit vector of (5 

value of singular integral 

kernel function 

unit vector at downwash point, 
normal to velocity vector 
unit vector at doublet point, 
normal to velocity vector 

direction cosines of n 

direction cosines of n 0 
pressure 

position vector of doublet 
from inertial frame origin 
position vector of downwash 
point from inertial frame origin 



q source or doublet strength 

Rt rotor tip radius 

R s rotor root radius 

r distance of downwash point 

along the span 

r 0 distance of doublet along the 
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upper limit of spanwise panel 
lower limit of spanwise panel 

distance of doublet along 
span at singular point time 
field time 

velocity of rotor system, 
parallel to x-axis, positive in 
negative x-direction 

velocity at downwash points 

4 - 

velocity component of V at 
the downwash point normal to 
the rotor leading edge 

velocity of doublet 
velocity of rotor system, 
parallel to z axis 
downwash velocity 
distance from pitch axis to 
downwash point 
Cartesian coordinates of 
downwash point 
Cartesian coordinates of 
doublet position 
twist angle at downwash point 
twist angle at doublet 
position 

angle of axis of rotation 
relative to z-axis 


V 0 /c 




Vc= 


1_ 

C 3t 

angular position of blade at 
time t 

angular position of blade at 
time t 

blade angle of attack 

blade angle relative to plane 

of rotation 

advance ratio 

air density 

time 


time at which integrand in Eq. 

(24) becomes singular 

velocity potential 

source acceleration potential 

doublet acceleration 

potential 

azimuth angle 

rotation speed of rotor 

vibration frequency of rotor 


BASIC FORMULATION 

The formulation of the aerodynamic 
equations is based on the linearized 
acceleration potential, approach. The fluid is 
considered perfect, with no separation and the 
formulation is based upon the assumption of 
small perturbations. The wake created by the 


lifting rotor is assumed to lie in the skewed 
helical path taken by the rotor blade. One 
reason for adopting the acceleration potential 
approach is that the pressure discontinuity 
occurs only on the surface of the blade and thus 
the boundary conditions need only be applied on 
the blade surface and not throughout the wake. 
The blade is treated as a very thin surface of 
discontinuity across which a pressure jump 
occurs. The effect of compressibiity is taken 
into account by utilizing the complete 
linearized potential for a lifting doublet, 
along with the effects of retarded time. 

As shown in Fig. 1, an .inertial coordinate 
system has been used in which the origin of 
coordinates is fixed to a point on the ground. 
The helicopter rotor is moving in the negative 
x-direction with velocity U, in the positive 
z-direction with velocity W and is rotating 
counter clockwise with a constant angular 
velocity fl. A point of interest on the rotor 

4 - 

blade is designated by the radius vector X q (t) 

from the origin of the ground based coordinate 
system. 

Let Y be the acceleration potential of a 
source (or doublet), the perturbation pressure 
is then’ given by 


p = -pY (1) 

This expression represents the pressure p at 

4 - 

point X due to a single source 

4 - 

(or doublet) located at X 0 . The potential Y 
contains a constant "q“ which represents the 
strength of the source and thus the magnitude of 
the pressure. In this form, there is no 
boundary condition available to determine the 
value of “q" and the resulting pressure. 

Recourse can be made to the velocity potential, 
since the spatial derivative of a velocity 
potential represents a velocity. The 
relationship between the pressure and velocity 
potential for an inertial coordinate system is 

P = -Pjt (2) 

where 4? is the substantial derivative, 
dt 

Dropping out the second order terms 

and integrating with respect to field time 

results in 

t 

<t»(t) = / Y(t‘) dt 1 (3) 


The acceleration potential Y s satisfies 
the wave equation 

2 

/v - iy —J- = -4n f(X,t) (4) 

s c at 

4 - 

where f(X,t) is a source distribution. 
Furthermore, if an isolated source is 

4 - -4 4 - 

moving with velocity V 0 , then f(X,t) = 6(X - 

-V 

V 0 t) where 6 is the delta function. 



Using the Green's function formulation the 
acceleration potential expression for a moving 
source, ¥ s can be written as (Morse and 
Feshbach, 1978, p. 841) 
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4 

where X 0 (t) designates the position of the 

4 

source at time x, X is the position of the 

4 - 

field point at the time t, V 0 (x) is the 
velocity of the source point at time t, c is 
the speed of sound and q is the strength of 
the source. An auxiliary equation which relates 
the time interval (t - x) to the distance 
between the two points is 


required. This derivative is taken normal to 
the flight path at the location of the downwash 
point, as follows 
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which is usually referred to as the causality 
condition. Eq. (5) expresses the potential as 
an explicit function of x, and only through Eq. 


( 6 ) as an implicit function of t ans t. From 
Eq. (3), the velocity potential due to a moving 
source is 


t t 

4 > s (t) = / V (t 1 ) dt‘ = / 

— 00 - 00 


q(T' ) 

dt 1 

4*(D - 8 • $J t , {7} 


Eq. (9) gives the downwash at a field point 
(x, y, z, t) due to a doublet placed at a point 
( x o»yo» z o» T ) having a strength q. In 
order to represent a lifting surface such as a 
rotor, it is necessary to distribute the 
doublets over the lifting surface and integrate 
over the surface to obtain the downwash at a 
field point. If the downwash is known, the 
quantity "q" can be determined. Letting K be 
the expression on the RHS of Eq. (9), the final 
equation is 

w n = // K dA' (10) 


x 



where ff = ? - jt Q , D = |5j 

and dt‘ = [1 - dx' 

The quantities x‘ , t' and t, x satisfy Eq. 

( 6 ). 


By definition, the doublet velocity 

potential $r of a doublet aligned along n 0 
can be written as 


*D (t) = Tn^ *s (t) 


V * 0 *s = ' n o 
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Note that for incompressible flow, c + «, the 
first term 0 and the integral remains 
unchanged except for the upper limit where x = 
t. 


To obtain the final equation for downwash 
Aw n , a second directional derivative is 


where A' is the area of the rotor surface. 

The LHS, w n , represents the known boundary 
condition and is the velocity normal to the 
velocity vector at the downwash point. By using 
the no flow condition for the velocity 
perpendicular to the blade surface, the velocity 

component in the n direction is V n tane w or 

w n (r,t) = V n tane w = ]j K dA' (11) 

where V n is the velocity component of V at the 
downwash point and is normal to the rotor 
leading edge and ^ is the angle of attack. 

Thus the problem requires setting up a method of 
solution of Eq. (11) from which a value of q, 
the unknown doublet strength, can be determined 
which satisfy the known velocity boundary 
conditions w n . 

This represents a rather formidable 
computing task and the history of lifting 
surface theory even for non-rotatiny wings has 
centered on devising approximate methods to 
accomplish the integration in an economical 
manner. One method, termed the vortex lattice 
method, has been very successfully applied to 
aircraft wings, and is probably the more 
economical procedure of the many variants. This 
method was first demonstrated for the unsteady 
case by Runyan and Woolston (1957) and was later 
expanded by Albano and Rodden (1969). This is 
the method adopted in this paper and the 
application will be discussed later. 
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Specification of Coordinate System 


The blade has the cord C and length Rt - 
R s , R s being distance to the the root of 
the blade, Rt is the distance to the tip of 
the blade. Let the blade momentarily coincide 
with the coordinate system along the positive 
x-axis at t = 0 and execute a counterclockwise 
rotation with angular velocity fl while moving 
with velocity U along the negative x direction 
and velocity W along the positive z 

direction. Since the vortex lattice method has 
been adopted, the doublet point lies C/4 ahead 
and the downwash point lies C/4 aft of the 
section midchord. The position of the doublet 
point as well as the downwash point can be 
established as follows. The Cartesian 
components of the doublet position are 

x 0 = -Ux + r 0 cos(flx) - (C/4)sin(flx)cos oq 
yo = r 0 sin(flx)+ (C/4) cos(nx)cos oq 
z 0 = Wx + (C/4) sin cxq (12) 

where r 0 is the radial distance of the doublet 
along the span. With the substitution of 
C + -C, r 0 r, x + t the position of the 
downwash point is given by 

x = -Ut + r cos(nt)+ (C/4) sin(nt)cos a 
y = r sin(flt) - (C/4) cos(nt) cos a (13) 

z = Wt - C/4 sin a 

In Eqs. (12) and (13), the angles a,OQ are the 

twist angles of the velocity vectors 'f and $ 0 , 
respectively, defined by 

- W 

tan ° " U si n(flt)+ rfl 

W (14) 

tan “o = u sin(flx)+ r Q fl 

The reference plane defined by the doublets and 
downwash points is a twisted surface. From 

Eq. (12) the doublet velocity can be computed, 
namely the time derivative of the position 
vectors. 


V o = V + V + z o k 

The unit vector n 0 is chosen to be 
perpendicular to the twisted surface created 
-► 

by the velocity vector V 0 which is a function 
of r 0 , through Eq. (14) . 

-*■ 

Express n 0 as 

% = V + V + n o k (15) 

where l , m , n are the directional cosines of 
o o o 


*o = 


W(U + r Q n si n ( £ 2 x)+ (C/4)n cos(flx)cos a Q ) 


L /w 2 + v ' 2 


W(r Q fl cos (nx) - (C/4)n sin(nx)cos a Q ) 


m Q = - 


v:/w^7 


'(16) 


o j 

/w 2 + v ' 2 


where 

V o l2 = (U + r Q n sin(flx)+ (C/4)n cos(nx)cos o Q ) 2 
+ (r Q s cos(Kx) - (C/4) a sin(flx)cos a Q ) 2 (17) 


By the same procedure, n = Jti + mj + nk, 
where 


_ W(U + ra sin(nt)- ( C/4 ) a cos(nt)cos a) 

vV W 2 + v ' 2 

3 -W(rfl cos(nt)-t- (C/4)« sin(flt)cos a) 


(18) 


•A 2 + 


7L 


n = 


fFTT 1 


and 

V' 2 = (U + rs sin(nt) - ( C/4) a cos(nt)cos o) 2 
+ (rfl cos(flt)+ ( C/4) fl sin(flt)cos a) 2 (19) 

•> 4 - 

the vector 0 = X-X 0 defined in Eq. (7) can be 
expressed as 

D = {[U(t-x) + r cos(flt)- r Q cos(flx) 

o 

+ (C/4)(sin(flt)cos a+ sin(flx)cos a Q )] 

+ [r sin(flt)- r Q sin(nx) (2°) 

2 

-(C/4)(cos(nt)cos a + cos(nx)cos a Q )] 

+ [W(t-x) -(C/4) (sin a Q + sin a )] 2 } 1/2 

With the substitution of the quantities, the 
integral Eq. (11) was solved for the unknown 
q(r 0 >x) by using a collocation process based 


the unit vector rio . It can be shown that 
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on the vortex lattice assumption. The kernel is 
singular when 0=0, and this was handled by use 
of the finite part technique. 


SOLUTION OF INTEGRAL EQUATION 

In following the vortex lattice technique 
the rotor is divided into a number of 
predetermined panels, both spanwise and 
chordwise, In each chordwise panel, a line of 
doublets of unknown strength q^ is located at 
the 25% chordwise location of the particular 
panel, and the downwash is evaluated at the 
point located at 75% chordwise location of the 
panel. Therefore, a collocation procedure is 
used to obtain a set of equations in terms of 
the unknown loadings qj . It is also assumed 
that the spanwise loading qi is constant along 
each of the panels. A set of equations is thus 
obtained as shown below. 

w = l A q (21) 

n " nm^m 
m 

r 

where A = / u K dr and where n refers to 

nm 1 nm o 


singularity. The integration domain was divided 
into areas as shown in Fig. 2. Areas 1-4 
(hatched) were computed numerically using a 
two-dimensional Romberg integration (Davis and 
Rabinowitz, 1967) and the contribution of the 
singular region (unhatched) was obtained in 
closed form by consideration of the finite part 
as shown in the next section. 

Treatment of Singular Term in Integral - The 
integral in the downwash equation, Eq. (11), is 
singular when D+0 and produces a complication 
which must be treated properly. It should be 
remembered that the integration path along "x" 
is the path the doublet has taken in arriving at 
the final doublet point at (c/4, r 0 ) measured 
in the local blade coordinates and can be 
considered as the wake. The integration takes 
place along the path from -~ to the final 
doublet position at x 0 . The distance D is the 
distance from the integration point at time x to 

the downwash point at t. 

There is a particular set of values of r 0 
and x for which the denominator D approaches 
zero, thus resulting in an infinite integrand. 
The singular part of the Eq. (11) is 


the downwash point and m refers to the vortex 
lattice. The kernel K is a complicated function 
which involves an integration over x. 



-»■ + 
n* n 0 - 


3(D *n) (D • n Q ) 


dx dr Q (24) 


The term q(r 0 ,x) represents the strength 
of the doublet located at r 0 and at time x, 
and is proportional to the unknown loading. In 
order to account for unsteadiness, a solution 
was formulated to take into account the time 
variation of the strength of the wake. This was 
done by assuming a Fourier series of the form 
m 

q(r Q ,x) = A Q + £ (A n cos(nlix) + B n sin(nflx)) (22) 

If q(r 0 ,x) is assumed to be a function of r 0 
alone, which means that the wake strength does 
not vary with time, the Fourier series reduces 
to q(r 0 ) = A 0 . A solution obtained with 
this approximation is termed the quasi- 
steady solution. 

This series was inserted in the basic 
equation and integrated with respect to x. 
However, there were more unknowns than 
simultaneous equations to solve for the 
unknowns. The additional required equations 
were obtained by evaluating Eq. (11) at a number 
of azimuth locations. For instance if m = 1, 
then 

q(r ,x) = A„ + A, cos«x + B.sinflt. (23) 

The azimuth was divided into equal segments of 
120° and the proper boundary conditions 
applied at i|i = 0°, 120°, and 240° thus 
prQV.id.ina. the necessary additional equations. 

Numerical Integration of Kernel 

the integration was performed by numerical 
integration, except for the area surrounding the 
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As ft-H) at the downwash point, D becomes 

perpendicular to n, therefore, at the singular 
point, the second term is zero and will be 
neglected in the treatment of the singularity. 
However, this second term is retained in all of 
the numerical integrations involving Areas 1-4 
since it represents an important contribution 
particularly when the blade is passing over a 
trailing wake. 

The time and distance at which the integral I 

becomes singular are designated by x and r 0 . 

The domain of the integration in Eq. (24) 
consists of a rectangle in which the duration 
x 2 -xi is kept extremely small. In other 
words, the integration is performed along a slit 
in r 0 , over which the 2nd term in Eq. (24) is 
negligible. Therefore the integral I can be 
approximated by 

r u n • n 

I = / / 2 1 1 dx dr (25) 

T 1 D 

Furthermore, noticing that D 2 is quadratic in 
fo* °*o is independent of r 0 , then the 
integration on r 0 can be performed 
analytically. This can be achieved by 
recognizing that in the vortex lattice method, 
the rotor is divided into spanwise panels from 
r^ to r u . If these spanwise panels are 
small then the variation in oq is small. 

da o „ Wn 

dr o (U sine Q + r Q n) 2 + W 2 ^ 

If the value of oq is approximated by its mid 
panel value, it is possible to integrate Eq. 
(24), in closed form in the r 0 direction. 

This is quite acceptable in the helicopter mode. 



because dao/dr 0 is in the order of 

magnitude 10"3 or smaller. The value n 0 is 
also a function of r 0 , but in the region of 
the singularity it has a very small variation 
and Is evaluated at the singular position, 

A A 

(x,r 0 ). Performing the r 0 integration 
results in the form 

1 ■ $$ d ' <Z7) 


where g(x) is a function containing all the 
non-singular part after performing the r 0 
integration and f(x) = 0, at 

A A 

x=x ( ti < x < T 2 ) . It can be argued 
physically that since the quantity D(x,r 0 ; 
t,r) as well as its modified form f(x) (after 
integration over r 0 ) represents the distance 
between two points in space it must be positive 
and real for all its arguments, and never become 
negative. Denote the value of r 0 and t at 

which D becomes zero as r 0 and t 0 . Thus, in 

the neighborhood of x the function f(x) behaves 
like a parabolic function and has a second order 
zero. 

Expanding f(x) in a Taylor series about the 
singular point x results in 

f(x) = f(x) + f 1 ( x) ( x- x) + f"(x) (x - x) 2 /2 +... 

(28) 

Since x is a second order zero 

and f(x) = f’(x) =0 (29) 


Eq. (29) has been verified numerically, 
the square term is kept in Eq. (28), Eq. 
can be written as 

T ° L_ | 9(t) + £^1 + g"( x) . 


I = K 


1 f"(t) (x-x) 


(t-t) 


If only 
(27) 


dx (30) 


as Ax is kept large because the very large 
values of the integrand near the singularity are 
avoided. On the other hand, regarding the 
finite part integration, the denominator was 
expanded in a Taylor series about the 

singular point, x. Therefore, it is desirable 
to maintain Ax as small as‘possible to keep 
within the limits of the applicability of the 
series expansion. Numerous calculations were 
made, varying Ax until a reasonable 
convergence was found. This value was found to 

be .01 (t- x) , i.e. 1% of the time difference. 
Actually, there is very little difference 
between 1% or 10% of the time difference and the 
computing time and cost is considerably reduced 
by using 10%. For trend studies 10% is 
recommended principally to reduce computer 
costs, However, for final design type analysis, 
a smaller value of time difference Ax is more 
appropriate. 

For the spanwise direction, Ar 0 is also 
an integration limit variable. The finite part 
integral was obtained by approximating the angle 
of twist of the velocity vector across a segment 
by assuming it constant across the segment, 
having 'a value as determined at the center of 
segment. Numerical experimentation indicates 
that for a helicopter, Ar 0 = 0 is satis- 
factory. 


APPLICATION TO SPECIFIC EXAMPLES 

The foregoing analysis has been applied to 
several specific examples which are given in 
Figs. (3) and (4). The following section 
presents results for several paneling 
configurations; e.g. 5 spanwise and 1 chordwise 
panels (designated (5-1)) and 7 spanwise and 3 
chordwise (designated (7-3)). The rotor blade 
was maintained at a constant pitch setting of 
6fj = .1 radians for all the calculations. 


In Eq. (30), if x 2 and xj are chosen 


Single Blade 


symmetrically about x, then the odd derivative 
terms integrate to zero. Futhermore, the third 

term can be neglected since g"(x) is small. The 
major contribution comes from the. first term. 
Then using the standard integration technique 
(Mangier, 1952) the final result for the 
integral is 


I = . Mil 

f"(T) AT 


(31) 


where 2Ax ■ xg - xj and x^ < x < X 2 . 

A numerical problem arises because the 
finite part integration results in a negative 
number which is close to the total of the 
surrounding numerical integration areas which 
are positive. Thus, it is necessary to take the 
difference between large numbers, and the final 
integration accuracy is dependent on the 
accuracy of the two integrations. On the one 
hand, the numerical integration is more accurate 
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In order to investigate the convergence of 
the method when using the vortex lattice 
procedure, the program was run for several 
chordwise and spanwise elements for the incom- 
pressible case. The thrust coefficient Ct 
vs. the azimuth angle is shown in fig. (5). (In 
all of the following plots for thrust coeffi- 
cient vs. azimuth angle, the thrust was 
calculated for 16 uniformly spaced azimuth 
angles and each curve was faired using a 
cubic spline). The rotor was first divided into 
5 spanwise and one chordwise (5-1) panel and the 
results are shown by the solid line. The 
chordwise division was increased to (5-2) and 
the results are shown by the long dashed line. 

It can be seen that very little change has taken 
place. The spanwise divisions were increased to 
(7-1) and the largest change occurred at 
t = 0° where the difference in Cj is about 
11%. Increasing the chordwise divisions to 3 
(7-3) shows convergence of the (7-1) case to be 
very good . 

An interesting phenomena occurs in the 
region of small azimuth angles. For ip=0 to 


* 


ft 


V 



37°, the lift increases to a local maximum at 
i(i= 370 then the lift abruptly falls to a local 
minimum for' y=60° and then rapidly increased 
to a maximum at ^=100°.’. A similar phenomenon 
is shown analytically by Eglof and Landgrebe 
(1983) in Fig. 60 of that report where a 
local minimum and a local maximum occur in the 
same ranye of azimuth angles, even though th.e 
geometry of the two blades and the flight 
conditions are different. Also, in Fig. 93 of 
the same report some test data shows a similar 
variation of loading in the same azimuth ranye. 

The chordwise pressure distributions for 
the (7-3) case are presented in figure 6. It 
should be remembered that in using the vortex 
lattice method, the loading is concentrated at 
the location of the vortex which for the (7-3) 
case is located at .0833C, .416C, and .75C. The 
pressure was faired using a cubic spline through 
the three vortex locations and the known value 
of zero at the trailing edge. The distributions 
are given for 7 spanwise positions. In general, 
the curves exhibit the expected shape, having 
the largest values as the leading edge is 
approached. For the span distribution the 
values at r/Rj = .85 are slightly larger than 
the values at r/RT = .95,, indicating a falling 
off in the tip region. 

From these concentrated forces, the section 
pitching moment can be calculated. Figure 7 
presents these results for f = 90 degrees. The 
section moment was taken about the 1/4 C and a 
nose down moment is taken as positive. The 
pitching moment shows some rather dramatic 
changes along the span. The moment is nose up 
near the tip (r/Rj = .95), changes to a small 
nose down value, then becomes nose up for most 
of the inboard reyion. Integration of the 
moment would result in a total pitch moment up 
at f * 90°. 


Swept Tip 

The segments used for the vortex lattice 
for the swept tip studies were(5-l), where two 
equal segments were used in the tip region and 
three equal segments were used in the unswept 
inboard section. In Fig. 8 the lift is shown 
plotted against azimuth for the two sweep 
conditions and for zero sweep. In general, the 
three results show little difference. The 
sweptback configuration has a larger lift from 
t » 300° to 40°. For y « 100° to 240°, 
the swept forward configuration has a sliyhtly 
larger lift. It appears that the total lift for 
one rotation for the swept-back case and the 
sweptforward case would give about the same lift 
as produced by the unswept rotor. In Fig. 9 
the lift distribution along the rotor span is 
given for 1 ) 1 = 0°. The major effect of sweep 
is concentrated at the tip, where the swept-back 
tip load is greater than both the unswept and 
sweptback cases . In fig. 10, tf ,, = 180°. 
Comparing to fig. 9, the swept-back tip load is 
larger than both the unswept and the 
swept-forward tips. 


Blade Oscillating in Pitch 

An example of unsteady loads oh a rotor 
blade with (5-1) paneling which is oscillating 
in a pitching mode about the mid-chord at a 
frequency of 4 per revolution (120 cycles/sec) 
is given on fig. 11. For this case a 17 term 
Fourier series (m=8) was used to simulate the 
oscillating load, which was comprised of one 
constant term, 8 cosine terms; and 8 sine 
terms. The steady and unsteady rotor blade 
loading is yiven for one revolution. The blade 
was oscillated through an angle of .1 rad. 
about a mean angle of .1 rad. The effect of the 
oscillation is readily apparent as compared to 
the steady case. With the harmonic 
representation of the loading, the magnitude and 
phase of the several harmonic loads are easily 
determined. The magnitudes are plotted in 
Fig. 12. The only harmonic loads that were 
significantly changed from the steady case were 
the 3rd, 4th and 5th. Both the 3rd and 5th 
harmonics were increased and the 4th harmonic 
was dramatically increased. Another calculation 
was made for the non-oscillatory unsteady case 
and compared to the quasi -steady case. 

Virtually no difference was observed, indicating 
that, at least for this case, the rate of change 
of loading in a revolution of the blade is small 
enough so that the effect of a variable wake is 
negligible. 


Compressible Effects (5-1) 

For a one-bladed rotor, the effect of 
compressibility is illustrated in Fig. 13, in 
which the Cy is plotted ayainst azimuth 
angle. The incompressible result is included 
for comparison. As expected, the compressible 
load is larger than the incompressible 
throughout one revolution. The effect is 
yreatest in the reyion of the advancing blade 
and smallest in the retreating reyion as would 
be expected. 


Two-Bladed Rotor in Compressible Flow (5-1 per 
bl ade) 

The method has been extended to the 
two-bladed rotor for the compressible case and 
the results are shown in Fig. 14. The thrust 
coefficient Cj per blade is yiven vs. azimuth 
angle for a single bladed rotor and for a 
two-bladed rotor. For azimuth angles from 
ip' = 20° to 120° the single blade rotor has a 
larger Cj. For <l> = 120° to 260°, the Cq- 
on the one and two-bladed rotors are 
approximately the same. However, for i|j = 260° 
to 340° a dramatic reduction in lift occurs 
for the two-bladed rotor as compared to the one 
bladed results. The lowest lift occurs at 
<li = 292° which places the other blade of the 
two-bladed rotor at 112°, the point of 
maximum lift on the other blade. Apparently the 
hiyh lift on the blade at <l>= 112° creates a 
very unfavorable induced velocity on the second 
blade at 0= 292° which requires the loading 
to yo to zero in order to satisfy the boundary 
conditions at <p = 292°. 
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CONCLUDING REMARKS 

A linearized lifting surface theory 
including the effects of compressibility has 
been developed for a helicopter rotor in forward 
flight. The method utilizes the concept of the 
acceleration potential, and makes use of the 
vortex-lattice procedure for performing the 
required integrations. In addition, the method 
has been extended to include the effects of 
unsteady flow. 

Sample calculations have been done for 
several cases. These include the effect of 
swept-back and swept-forward tip. The effect of 
these two tip configurations was minimal on the 
total loading for one revolution. However, the 
loading distribution changed considerably for 
several azimuth positions. A comparison of the 
thrust coefficient, Cj, of a one bladed rotor 
and a two bladed rotor was made. In the 
azimuthal range between 20° and 120°, the 
one bladed rotor showed higher lift. However 
between = 260° to 340° the two bladed 
rotor indicated a lower Cj. Compressibility 
was investigated for one configuration. As 
expected, the effect was greatest in the 
advancing blade region ( 90°) and was 

minimal in the retreating blade region. The 
effect on Cj of a blade oscillating in pitch 
at 4/rev is given. The effect on the total 
blade lift is shown and the effect of the 
oscillation is readily apparent. The harmonic 
content was calculated and the greatest 
difference between the oscillatory and 
non-oscillatory cases was found in the 4th 
harmonic . 
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Fig. 1 Inertial Coordinate System 



Fig. 2 Integration Areas 
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Fig. 3 Unswept Configuration and Input 
Parameters 


Fig. 4 Swept Tip Configurations 




Fig. 5 Thrust Coefficient vs. Azimuth Angle for 
Four Panel Configurations for a Single 
Rotor Blade, incompressible, (u = 0.17, 
00 = 0.1 rad, a r ■ .OS rad, 
fl = 30 rad/ sec) 


Fig. 6 Chordwise Pressure Distribution for 
Several Spanwise Locations, ip = 90°, 
incompressible, (u = .17, eg = 0.1 
rad, o r = .05 rad, a = 30 rad/ sec) 



Fig. 7 Spanwise Section Moment Distribution 
about 1/4C - Positive Nose Down, 

= 90°, incompressible, 

(u = 0.17, eg = 0.1 rad, o r = .05 
rad, a = 30 rad/ sec) 


Fig. 8 Comparison of Lift on a Swept-Back Zero 
Sweep and Swept-Forward Blade, 
incompressible, (u = 0.17, e 0 = 0.1 
rad, o r = .05 rad, n = 30 rad/ sec) 
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Fig. 9 Spanwise Section Lift Distribution for 
Swept-Tip Configurations, f = 0°, 
Incompressible, (u = 0.17, e 8 = 0.1 
rad, a r = .05 rad, ft - 30 rad/sec) 


• Fig. 10 Spanwise Section Lift Distribution for 
Swept-Tip Configurations, <|> = 180°, 
incompressible, (y = 0.17, eg = 0.1 
rad, a,. = .05 rad, ft = 30 rad/ sec) 



Angle, deg 

Fig. 11 Comparison of Lift on a Rotor Blade 

Oscillating in Pitch at 4/Rev. to the 
Lift on a Non-Oscillatory Blade, 
incompressible (e =0.17, eg = 0.1 
rad, a r = .05 rad, a = 30 rad/sec) 



Fig. 12 Harmonic Content for Non-Oscillatory and 
Oscillatory Cases - r/Ry = .95, 
incompressible, (y = 0.17, 6g = 0.1 
rad, a r = .05 rad, ft = 30 raa/sec) 
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Fly. 13 Incompressible and Compressible Lift for Fig. 14 Lift on Two-Bladed and One-Bladed Rotor 
One-Bladed Rotor, (p = 0.17, 6a = 0.1 vs. Azimuth Angle, compressible, 

rad, a r = .05 rad, Mjjp = 0.54) ( p = 0.17, 6o = 0.1 rad, « r = .05 

rad, Mjjp = 0.54) 
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